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CsJ I Abstract 

In this paper we studied the asymptotic eigenvalue statistics of the 2 matrix model 
with the probability measure 

Z,; 1 exp (-n (tr(V(Mi) + W{M 2 ) - tM 1 M 2 )) dMidM 2 , 

'■^ ■ in the case where W = ^ and V is a general even polynomial. We studied the 

correlation kernel for the eigenvalues of the matrix M\ in the limit as n — > do. 
We extended the results of Duits and Kuijlaars in [14J to the case when the limiting 
eigenvalue density for M\ is supported on multiple intervals. The results are achieved 
by constructing the parametrix to a Riemann-Hilbert problem obtained in [T3] with 
theta functions and then showing that this parametrix is well-defined for all n by 
studying the theta divisor. 

o 

1 Introduction 

q . 

1.1 2 matrix models and biorthogonal polynomials 

oo : 

O ■ The 2-matrix Hermitian models are matrix models with the probability measure 

; Z~ l exp (-n {tr(V(M 1 ) + W(M 2 ) - tM 1 M 2 )) dMxdM 2 , (LI) 

H 

denned on the space of pairs (M 1; M 2 ) of n x n Hermitian matrix. The constant Z n is the 
normalization constant of the measure, rGl \ {0} is a coupling constant and dM 1; dM 2 
are the standard Lebesgue measures on the space of Hermitian matrices. In (11. ip . V and 
W are called potentials of the matrix model. In this paper, we shall consider V to be a 

4 

general even polynomial and W to be the monomial W(y) = \. 

Let X\, . . . ,x n and yi, . . . , y n be the eigenvalues of the matrices M 1 and M 2 respectively, 
then the eigenvalues of the matrix model (11 .11) are distributed according to 

= Zn 1 fl(x, - x k )\y t - y k f e ~<^Un^M-^) (1.2) 

i<k 
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where Z n is a normalization constant and x — (x\, . . . , x n ), y = (yx, . . . , y n ). 

The two-matrix model was introduced in [27] , [31] as a generalization of the one matrix 
model to study critical phenomena in physical systems. The 2 matrix model is needed 
to represent all conformal models in statistical physics [9]. It is also a powerful tool in 
the studies of random surfaces as the large N expansion of the partition function Z n is 
expected to be the generating function of discretized surface [29]. Since its introduction, 
the 2 matrix model has become a very active research area p], [2], [3], [I], [5], [6], [13], [33], 
P21, PS], [20], [19], [21], [22], [23], [28], [30] and one of the major problems is to obtain 
rigorous asymptotics for the eigenvalue statistics. A good review of the subject can be 
found in [T5] , [To] . 

A particular important object in the studies of eigenvalue statistics is the correlation 
function tl^Xi, . . . , x m , y u . . . , y,) 

/ j\2 f f n n 

Ku(* 1; . . . , x m , Vl , . . . , Vl ) = (w _ m "i (w _ , J mn _ m V(x, y) jg be, ft dy k . 

(1.3) 

In [18], [32], a connection between biorthogonal polynomials and the correlation func- 
tions of 2 matrix models (11. 3ft was found. Let pk{x) and qi(y) be monic polynomials of 
degrees k and / such that 

dxdy Pk (x) qi (y)e- n ^ +w ^- T ^ = h k S kh (1.4) 

for some constants h k , then these polynomials exist and are unique [3], [TTJ. These poly- 
nomials are known as biorthogonal polynomials. 

Let us define some integral transforms of the biorthogonal polynomials by 



;i.5) 




Q k (x) = e- nV ^ / q k (y)e- n ^- T ^dy, 
Pk(y) = e~ nW{y) [ Pk (x)e- n ^^dx, 



and define the kernels to be 



n— 1 . n— 1 ^ 

^u(^i,x 2 ) = ^—p j (x 1 )Q j (x 2 ), K% 2 (x,y) = ^2—P j (x)q j (y), 
n— 1 

3=0 3 
n—l 

3=0 3 
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Then the correlation function l (11.31) has the following determinantal expression. 

Upon averaging over the variables we see that the m-point correlation function TZ^ nQ {x\, . . . , x m ) 
for the eigenvalues of the matrix Mi is given by the kernel K^Xi, Xj), 

n n mfi (x u ...,x m ) = det (K^( Xl ,x 3 ))^ =1 (1.8) 

The purpose of this paper is to provide a rigorous asymptotic expression for the kernel K™ x 

4 

as n — ► oo for W(y) = v and V(rr) being a general even polynomial. 

Due to a generalized Christoffel-Darboux formula, the kernels K™ x can be expressed in 
terms of a finite sum of the biorthogonal polynomials (See [5], (TH], [33] and [7]). This re- 
duces the problem of finding an asymptotic expression for into finding the asymptotics 
of the biorthogonal polynomials. 



1.2 Rigorous results in the "one-cut regular" case 

Until recently, most results in the asymptotics of biorthogonal polynomials have been 
obtained through heuristic argument (See [TU], [20] )• For a long time, the only rigorous 
result was the case when both W(y) and V(x) are quadratic polynomials [IT]. In the recent 
work by Duits and Kuijlaars [H], (See also p3], Chapter 5, which was later made into the 
publication [Ti]). the Deift-Zhou steepest decent method ([TU]. [TT] . [T2] . see also [Ej) was 

4 

successfully applied to obtain the asymptotics of biorthogonal polynomials with W(y) = \ 
and V(x) an even polynomial in the case when the limiting eigenvalue density for Mi is 
supported on a single interval. The main idea in [14] is to transform and approximate the 
Riemann-Hilbert problem satisfied by the biorthogonal polynomials [30], [6] (See Section [3] 
for details) via the use of suitable equilibrium measures and then solve the approximated 
Riemann-Hilbert problem explicitly to obtain asymptotic formula for the biorthogonal 
polynomials. The results in [H] was obtained in the case when one of these equilibrium 
measures is supported on a single interval. This case is called the "one-cut regular case" 
in P3|. 

To be precise, let I(ut, Vj) be the following energy function 

Hui,Vj) = j Jlogfj^^Adu^dUjiy), (1.9) 

where the integral is performed on the supports of the measures ^ and Vj. Then the 
equilibrium measures /ii, \ii and /X3 are the measures that minimize the following energy 
function Ey(vi, u 2 , vz). 
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Definition 1. (Definition 2.2 in [Lflj (See also Definition 5.2.1 in fTB^)) The equilibrium 
measure fi 2 , ^3) is the triplet of measures that minimizes the following energy function. 

E v {v 1 ,v 2 ,v z )=^I{v h v j )-^I{v h v j+1 ) + I (v(.)-^|z|iW). (1.10) 

amount non-negative Borel measures V\, v 2 and U3 that satisfy the following properties. 

1. All the measures Vj, j = 1, 2, 3 have finite logarithmic energies; 

2. v\ and z/3 are measures supported on M with i^i(R) = 1 and v% (R) = |; 

3. v 2 is a measure supported on iM with u 2 {iM) = |; 

4. Let a be the unbounded measure on iM. given by 

/3 

da(z) = y— T *\z\s\dz\, zeiM, (1.11) 
2n 

where \dz\ is the unit arc length on iM, then v 2 satisfies the constraint v 2 < a. 
Let U u (x) be the logarithmic potential of the measure v. 

U u {x) = - J \og\x - s\dv(s), (1.12) 

then it was shown in [T4] that the logarithmic potentials of \i\ and \i 2 satisfy the following 
properties 



2C/« ( x ) = U*» (x) - V (x) + -T3 \x\ 3 + I, x G S, 



Ml ' 

(1.13) 



o 

2U^(x) > U» 2 (x) - V{x) + - T s\ x \s + l, xeM\S^, 

for some constant /, where is the support of The properties of the equilibrium 
measures fi\, \i 2 and were studied in [TJ] and we have the following 

Theorem 1. (Theorem 2.3 in [14] (see also Theorem 5.2.2 in [13])) Let V be an even 
polynomial and r > 0. Then there is a unique minimizer (/ii, [i 2 , ^3) of E{i>i,v 2 ,Vz) in 
$1.101) that satisfies the conditions in Definition d Let us denote the support of the Borel 
measure v by S u , then we have 

1. consists of finitely many disjoint intervals 

=Uf+ 1 [A 2j _ 1 ,A 23 -], (1.14) 
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where Xj € R and the points are ordered such that Xj < Xk if j < k. Moreover, pi 
is absolutely continuous with respect to the Lebesgue measure, and on each interval 
[A^j-i, X2j], it has a continuous density of the form 



Pi{z) = pi(z)dz = ipj{z)J (X 2 j - z)(z - X 2j -i), z E [X 2 j-i, X 2 j], (1-15) 
where ipj(z) is analytic and non-negative on [X 2 j-i, X 2 j]. 

2. Let o be the measure in U.ll\) , then = iR and there exists a constant c > such 
that the support S a ^.^ 2 of a — p 2 is given by 

S a - tla =m.\(-ic,ic). (1.16) 

Moreover, a — p 2 has an analytic density on S ff - m that vanishes as a square root at 
He. 

3. = R and p% has a density which is analytic in R \ {0}. 

4. For j = 1,2,3, we have fij{A) = fij(-A) for any Borel set A. 

Remark 1. In particular, by 4- in the above, we see that S Ml is symmetric under the map 
z 1 — > — z, that is, we have Xk = —X 2g + 2 -k+i- 

We then have the following definition of regularity. (See Definition 2.5 in [H]) 

Definition 2. The potential V(x) is regular if the following conditions are satisfied. 

1. The inequality in $1.13\) is strict outside of S^; 

2. The density p\ (z) vanishes like a square-root at the end points of 5* Ml ; 

3. The density p\ (z) does not vanish in the interior of . 

It is known that a generic potential V(x) is regular p3]. The "one-cut regular case" is 
the case when consists only of a single interval and that V(x) is regular. In [Hj, rigorous 
asymptotics of biorthogonal polynomials was obtained for this case. The asymptotics of the 
biorthogonal polynomials were then used to obtain an asymptotic expression for the kernel 
i^n in (jl.6p . In this paper we will extend these result to the case when pi is supported on 
any number of intervals. (See Theorem [2] and Theorem [3]) 

2 Statement of results 

In this paper we obtain universality results for the 2 matrix model with potentials W(y) = 

4 

and V(x) a general even polynomial. Moreover, we will assume the potential V(x) 
satisfies the regularity condition in Definition [2] and that n is a multiple of 3. Then we 
have the following result on the global eigenvalue distribution of the matrix Mi. 
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Theorem 2. Let (/xi, //2, ^3) be the equilibrium measures that minimize the functional 
l\1.10\) . Then as n — > 00 and n = 0(mod3) ; we have 

lim -K? 1 (x,x) = p 1 (x), (2.1) 

n-^oo 

uniformly for x6K, where p% is the density of \i\ in Definition [H 

As explained in [T3j, the requirement n = 0(mod3) is not essential and is only imposed 
to minimize the technicality. The other result is the universality property of the kernel 

Theorem 3. Let (/ii, /i2, Ps) be the equilibrium measures that minimize the functional 
U.10\) . Then as n — > 00 and n = 0(mod3) ; we have the followings. 

1. Let x* be a point in the interior of the support S^. Then we have 

hm x* + — , x* + — = V 2 - 2 

n^oonpi(x*) \ npi(x*) np\(x*) J n{u — v) 

uniformly for u, v in compact subsets o/R. 

2. Let (fj > be such that 

px{z) = ^±\z - A,-|5 +0(z- Xj) , 

7T 

as z — > Xj, j = 1, . . . , 2g + 2 inside of S^, where Xj are defined as in fll.l^| ). Taen 
we aave iTie following 

lim - — -/Y-', I A, + (-1V'— ^ A,, (-] ! 



™(n^)3 \ (n^-) 5 W 5 / (2.3) 

Ai(n)Ai'(n) - Ai(n)Ai'(n) 

u — v 

uniformly for u, v in compact subsets o/R, where Ai is the Airy function. 

Recall that the Airy function is the unique solution to the differential equation v" = zv 
that has the following asymptotic behavior as z — > 00 in the sector — 7r + e < arg(z) < 71 — e, 
for any e > 0. 

1 2 3 / 3 \ 

M(z) = r e"3 22 l + 0(z~5) \ -7T + e < arg(z) < tt - e, 2^00. (2.4) 



20Fz 
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where the branch cut of z? in the above is chosen to be the negative real axis. 

Although the steepest decent analysis in [2] already covers the general case without the 
1-cut assumption, solution to a 'modeled Riemann-Hilbert problem' (See (13.71) ) must be 
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obtained to complete the Riemann-Hilbert analysis and to extend the universality results 
to the general case. The main difficulties are to show that a solution of the modeled 
Riemann-Hilbert problem exists for all n and to find an explicit expression of it. This 
involves the study of the theta divisor, which is the set of points in which a theta function 
vanishes. (See Section H] for a more detailed description of the theta function) This is a 
difficult problem with very few results available. In this paper we managed to construct 
the solution of the modeled Riemann-Hilbert problem (13.71) with the use of theta functions 
and by using results from [21] and [25], we were able to show the existence of the solution 
M(z) to (13. 7p for all n. This allows us to extend the universality results in [H] to the case 
when V(x) is a general even polynomial and obtain Theorem [2] and Theorem [HI 

In many applications of the Deift-Zhou steepest decent method, theta function is needed 
to solve a modeled Riemann-Hilbert problem and the solvability of these modeled Riemann- 
Hilbert problems is important to guarantee the validity of these asymptotic formula for all 
n, as n — > oo. We believe the techniques and results in this paper will not only be valuable 
to the random matrix community studying 2 matrix models, but it will also be important 
to many other problems in which the Deift-Zhou steepest decent method is applicable. 
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3 Riemann-Hilbert analysis 

In this section we will summarize the results in [H] that are relevant to our analysis. We 
will not repeat the analysis in [H] but will simply state the results that are applicable to 
our studies. 

3.1 Riemann-Hilbert problem 

The biorthogonal polynomials Pk{x) in (jl.4p satisfies a Riemann-Hilbert problem [B], [30] 
similar to the one that is satisfied by orthogonal polynomials [26J. This allows the imple- 
mentation of the Deift-Zhou steepest decent method. ([10], [H], P2], [H]) 
Let Wj(x) be the weights defined by 

Wj (x) = e- nV ^ [ y j e~ n V r ~' n Vdy, j = 0,1,2. (3.1) 
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Assuming n is divisible by 3 and consider the following Riemann-Hilbert problem 
1. Y{z) is analytic in C \ R, 



Yjz)= Yjz) 



3. Y(z) 



fl 


w { 


(*) 


wx(z) 


w 2 (z)\ 









1 


















1 





1 


\o 










1 1 








(z n 





\ 


z' 1 


)) 





n 

Z 3 











z-f 











z~ 


n 

V 



2 e 



(3.2) 



oo. 



This Riemann-Hilbert problem has a unique solution given by the biorthogonal polynomial 
Pk{x) and some other polynomials, together with their Cauchy transforms [30J. 



Y(z) 



( p n (z) C{p n w ){z) C{p n wi){x) C{p n w 2 ){z) \ 

P { nUz) C( P W 1W0 )(Z) C(pW lWl )(z) C(pM lW2 )(z) 



VnUz) C( P ^ IW0 )(Z) C(tii lWl )( Z ) C( P ^ 1W2 )(Z) 



(1) 



\p ( n\(z) C(ti> lWo )(z) C(ti> lWl )(z) C(J>^2){Z)J 



(2) 



,(2) 



(3.3) 



where p^li(z), j = 0, 1,2 are some polynomials of degree n — 1 and C(/) is the Cauchy 
transform of the function /. 



C(f)(x) 



1 

2m 



■ds. 



s — X 



In [2], the Deift-Zhou steepest decent method ([ID], [H], [12], [H]) was extended to the 
Riemann-Hilbert problem (13. 3p . With the help of the equilibrium measures fix, fi 2 and 
/X3 that minimize fll.lOp . the authors of [H] were able to transform and approximate the 
Riemann-Hilbert problem (13.31) by an explicitly solvable one. To state their results, let us 
assume V(x) is regular. Let us denote an interval in the support ( 11. 14ft of /ii by Ej 
and a gap between the intervals by Ej. 



[A2J-1, X 2 j], j = 1, 

[\ 2j , \ 2 j+i], j = 1, 

(-00, Ai], E g+1 = 

We will define ctj to be the constants 



■ — 'j 
So 



...,0+1, 
■■■,9, 

[X 2g+2 , 00) 



(3.4) 



^1 ( u fife+i s i) > 



fc = 0, 



a 



9+1 



0. 



(3.5) 



Note that, since V(x) is an even polynomial, by Theorem [T], we have fii(A) = fii(-A) for 
any Borel set A. Therefore the constants at in (13. 5ft satisfy the symmetry 



g+l-k- 



(3.6) 
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Let us define the following Riemann-Hilbert problem for a matrix M(z). (See Section 8 of 
[H] and 5.10 of [13]) 

1. MO) is analytic inC\(RU <S CT _ M2 ), 

2. M+(z) = M.(z)J M (z), zeRU Sa-n, 

(I \ 
0^00 A 
10 \0 A 

\o o o z~y 

uniformly as z — > oo in the j th quadrant, 



3. M(z) = (J + 0(z~ x )) 



(3.7) 



4. M(z)=o((z-A i )"4), 
M(z) = O ((^T«c)^) , 



A 



j = l,...,2</ + 2, 



where S a ^^ 2 is oriented upward, the branch of is chosen such that Z3 £ R for z e M+ 
and the branch cut is chosen to be the negative real axis. The Aj are given by (with 



2wi 

uj = e s j 




w 1 w 
Ao = — \ 1 1 1 



v/3 



v/3 



v Co> 2 1 

-1 u 2 

-1 1 

-1 u - 



The jump matrices Jm(z) in f)3.7p are given by the followings 



Jm(z) 



Jm(z) 



Jm(z) 



( 1 0\ 
-10 
1 

yo o-io/ 

/ e -2nnia k g 




\ 

/l 0\ 

0-10 

10 

\0 1/ 



1 



0\ 

i 

i 1 ' 
-1 o) 

z £ Sfj—^ . 



z e Hjfc, fc = o, ... ,0 + 1, 



(3.8) 



The steepest decent analysis in [Hj leads to the 'modeled Riemann-Hilbert problem' ( 13.71) . 
Provided a solution M(z) of (13.71) exists and is uniformly bounded in n away from the 
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singularities, the analysis in [2] that leads to the asymptotic forms (12. II) . (12. 2ft and ( 12. 3D 
can be carried out with the parametrix M(z). 

Theorem 4. Let e > be a fixed small number independent on n. Let B £ j and B £ ,±i C be 
small discs of radius e centered at Xj and ±ic respectively. Let K, C C be a compact subset 
in C and let T be the set 



T=K\ (ujffBej U B £ , c U B e -ic) . 



(3.9) 



Suppose the solution M(z) of the Riemann-Hilbert problem Jff. 7| ) and its inverse M l (z) 
exist and satisfy the following conditions. 

1. Both M(z) and M~ 1 (z) are bounded in n uniformly inside T for any compact subset 

2. For any r > max{c, \2g-h2} independent on n, there exist constants Cj fc and w l ^ k , 
1 < < 4 ; independent on n, such that, for z > \r\, 



(M(z)) jk \<C l jk 





1 




Z$ 



for z in the I th quadrant, 
, for z in the I th quadrant. 



Then as n — > 00 and n = 0mod3, the asymptotic behavior of the kernel K™ x is given by 
( fH|) . (E3) and / TP) . 

In the following sections we will construct the solution M(z) with the help of theta 
functions and we will show that the solution satisfies the conditions in Theorem H] in 
Section El 



4 Theta function and Riemann surface 

We will now construct a Riemann surface from the equilibrium measures and use the 
theta function on this Riemann surface to construct a parametrix for the Riemann-Hilbert 
problem (13.71) . 

The Riemann surface is realized as a four-sheeted covering of the Riemann sphere. 
Define four copies of the Riemann sphere by Lj , j = 1 , . . . , 4 

Ci = C\S ftl , £ 2 = C \ (S^ U Str-nz), 
C 3 = C\(S a ^ l2 US, 3 ), £ 4 = C\(S M3 ), 

where C is the Riemann sphere obtained by adding the point z = 00 to C. 

The Riemann surface C is constructed as follows: L\ is connected to £2 via S^, £2 is 
connected to C 3 via S , cr _ M2 and £ 3 is connected to £ 4 via S^, as shown in Figure [TJ Let us 
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Figure 1: The sheet structure of the Riemann surface C. 



define the functions Fj(z) by 



s z 

3 



-dfij(s) 



(4.2) 



Then we have the following result 



Lemma 1. (Lemma 5.1 in [Lflj , (see also Lemma 5.4-1 in fHHf)) The function £ : Uj =1 Cj 
C defined by 





z)+V'{z), 




e&; 




Fi{z) 


_ , , 4 1 
- F 2 (z) + T3Z3, 


z 


eC 2 , 


Rez > 0; 


Fi{z) 


- F 2 {z) - ri{-z% 




eC 2 , 


Rez < 0; 


< F 2 {z) 


- F 3 (z) - rs(-z)s, 




eC 3 , 


Rez > 0; 


M*) 


- F 3 (z) + 


z 


eC 3 , 


Rez < 0; 


Fs(z] 


4-zri 4 1 
+ e 3 T3Z3 ; 


z 


e £ 4 , 


Imz > 0; 


. F 3 (z] 


2iri 4 1 
+ e 3 T3Z3, 


z 


e C 4 , 


Imz < 0. 



(4.3) 



has an extension to a meromorphic function (also denoted by on C. The meromorphic 
function has a pole of order degV — 1 at infinity on the first sheet, and a simple pole at the 
other points at infinity. We shall denote the restriction of £(z) to the sheet Ck by £&(z)- 

The Riemann surface £ is of genus g. Let us define a set of canonical basis of cycle as 
in Figure El The figure should be understood as follows. The top left rectangle denotes 
the first sheet the top right rectangle denotes C 2 , the lower left one denotes C 3 and 
the lower right one denotes £4. A 6-cycle is a loop in Ci around the branch cuts that is 
symmetric with respect to the real axis, while an a-cycle aj consist of a path in the upper 
half plane in C\ that goes from S^+i to Ej (Ej is defined in (13.41) ). together with a path in 
the lower half plane in C 2 that goes from Ej to The loop formed by these 2 paths 

is an a-cycle. We will also choose these 2 paths such that their projection on the complex 
z-plane are mapped onto each other under complex conjugation. 
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Figure 2: The a and b cycle of Riemann surface C. 

We can now define the basis of holomorphic differential that is dual to this basis of 
cycle. 

Let dcuj be holomorphic differential 1-forms on C such that 



dujj = 5 jk , l<j,k<g. 



(4.4) 



a k 



The 1-forms dujj are known as the holomorphic 1-forms that are dual to the basis of cycles 
(a, b). 

Let the 6-period of these 1-forms be II j, 



duj = Uij, 



(4.5) 



bi 



then the g x g matrix II with entries is symmetric and Im(II) > 0. 



4.1 Theta function and its properties 

-> C associated to the Riemann surface £ and this choice of 



The theta function 9 : C 9 
basis is defined by 



0(s) := £ 



inn-Uri+2iirs-ri 



(4.6) 
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The theta function has the following quasi-periodic property, which will be important to 
the construction of the parametrix. 



Proposition 1. The theta function is quasi-periodic with the following properties: 
9(s + M) =0(s), 



9(s + UM) =exp 



2-/ | - ( \I.s) - ( M,^M 



where (•, •) denotes the usual inner product in C 9 . 

We will now define the Abel map on £. The Abel map u : C 



9(s), 



C 9 is defined by 



(4.7) 



u(x) = (u^x), . . .,U g (x)) J 



x 



Xf) 



r 



(4.8) 



where x is a point on C. We will choose x so that x is the point on L\ that projects to 
A2 S +2 in C We will denote this point by \\ g+2 - 

The composition of the theta function with the Abel map is then a multi- valued function 
from C to C. It is either identically zero or it has g zeros on C. The following lemma tells 
us where the zeros are. 

Lemma 2. Let D = J2i=i di be a non special divisor of degree g on C, then the multi-valued 
function 

9(u(x) — u(D) — K) 

has precisely g zeros located at the points d i} i = 1, . . . , g. The vector K = (Ki, . . . , K g ) T 
is the Riemann constant 



A 2 9 +2 



Recall that a divisor of degree m is a formal sum of m points (counting multiplicity) 
on the Riemann surface and that two divisors Di and D 2 are equivalent if and only if 
there exists a meromorphic function f(x) on C with poles at the points of D 1 and zeros at 
the points of D 2 . A divisor Yli=i di is special if there exists a non-constant meromorphic 
function on C with g poles at the points d±, . . . , d g . The condition of D being non special 
is equivalent to the condition that 9{u{x) — u(D) — K) does not vanish identically. 

Theorem 5. Let di,...,d g be g points on a Riemann surface and let the multiplicity of 
dj within these g-tuple of points be kj. Then the function 9(u(z) — J2 9 j=i u (d g ) — K) is 
identically zero if and only if there exist a function f(z) that has poles of order kj at dj for 
j = 1, . . . , g and holomorphic elsewhere. 
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This is a consequence of the Riemann-Roch theorem. In general, for a given g + I 
points (counting multiplicity) on a Riemann surface, there are I independent meromorphic 
functions with poles exactly at these points. This can be thought of as an extension of the 
Liouville's theorem. 

Let 4>(z) be the anti-holomorphic involution on L defined by 

<Kz):{z,t{z))^QS,£<d) ( 4 -9) 

where is the function on C given by Lemma [U 

Then by the definition of the cycles in Figure [2J we see that under the involution 0, we 
have 

<P(bj) = -bj, (f)(aj) ~ a j} j = 1, . . . , g, (4.10) 

where the symbol ~ means that 4>(aj) is homologic to aj. 

In particular, if we consider the holomorphic 1-forms dujj(<f>(x)) on £, we have 



(p duj(4>(x)) = (b du)j(x) = (p dujj(x) = 5jk- 

J a k J<t>{flk) J a k 

Hence by the uniqueness of holomorphic 1-forms that is dual to the cycles (a, b), we have 
dujj((p(x)) = u)j(x). By computing the 6-periods of dwj(<j)(x)) and making use of (14.10p . we 
obtain the following. 

Lemma 3. The period matrix U of £ is purely imaginary. 

In particular, by (14. 6p . we see that ^(0) is real and positive and from Lemma [21 we 
see that 6{u{x)) has g zeros on C Let us denote these zeros by l%, ■ ■ ■ , t g . Then by the 
following result in [24j, we can simplify the expression of the Riemann constant K. 

Proposition 2. (See p. 308-309 of \24\) Suppose 6(u(x)) is not identically zero. Then it 
has g zeros in C. Let Li, . . . , i g be its zeros, then the Riemann constant is given by 

K = -Y^u{ h ). (4.11) 

j=i 

Remark 2. Let 5^ be the intervals in L\ that projects onto the gaps Ej in \3.4\ ), then as we 
shall see in CorollaryUl there is exactly one point Lj that belongs to H] US^ for j = 1, . . . , g. 
We shall label the Lj such that Lj G U . 

We would like to express the function 6(u(x)) as a meromorphic function on C with 
jump discontinuities. To do so, we need to define the contour of integration in the Abel 
map ( 14. 8 p such that the integral can be defined without ambiguity. We will define the 
contour of integration as in Figure El 

For a point z in C\ (£2), the contour of integration Si goes from \2 g +2 to z in C\ (£2) 
without intersecting (—00, \2g+2) an d the branch cuts on the imaginary axis. For a point 
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Figure 3: The contour of integration for the Abel map u. 



z in the upper (lower) half plane in £ 3 , the contour of integration S 3 consists of 2 parts. 
The first part lies in £ 2 , goes from A 2g+ 2 to a point i( (— i£) on the branch cut in the 
imaginary axis without intersecting (—00, A 2g+ 2) in £2 and enter the branch cut from the 
left hand side of this point. The second part lies in the upper (lower) half plane in £3, goes 
from the right hand side of i( to the point z. For a point z in the upper (lower) half 

plane in £ 4 , the contour of integration consists of 3 parts. The first part lies in C 2 , goes 
from \2 9 +2 to a point — i( on the branch cut in the imaginary axis without intersecting 
(—00, A 2£ ,+2) in £2 and enter the branch cut from the left hand side of this point. The 
second part lies in the lower (upper) half plane in £3, goes from the right hand side of — iC, 
(iQ to the origin. The last part lies in the upper (lower) half plane in £4, goes from the 
origin to the point z. The choice of the point ±i( in the construction is immaterial as long 
as it lies on the branch cut on the imaginary axis. 

Let z$ be the point on Cj that projects to z in C and A be a (/-dimensional constant 
vector. We can now define four functions 8^(u(z) + A) on the complex z-plane by 

9 j (u(z) + A) = 6{u{z j ) + A). (4.12) 

These functions will have jump discontinuities in the complex z-plane. By using the peri- 
odicity of the theta function ( 14.71) . we can compute their jump discontinuities. 
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Proposition 3. Let A = (A\, . . . , A g ) T be a g dimensional vector. The functions 8^(u(z) + 
A) are analytic mC\(lU SV-/^ ) . On R U S a -^ 2 , they satisfy the following conditions 



nl 


( ( \ 

M z ) 


+ A) 




J) 


+ A), 


z e E j} .7 = 1,.. 


6 + 


(u(z) 


+ A) 


= 9 l _(u 




+ A)e { ~ 


-l) l 2m(u j {z)+A J + ^f-^ 


9 + 


(u(z) 


+ A) 


= d l _(u 




+ A), 


z e (e u E g+1 j , 


el 


[u{z) 


+ A) 


= 6%{u 




+ A), 


Z G S a — i_i2 ; 


el 


(u(z) 


+ A) 


= e A T {u 




+ A), 


zGK, 




[u{z] 


+ A) 


= e l _( u (z) 


+ A), 


z G Sfj—^'i 1 1 



z e c 



1 = 1,2, 



j = l,...,g, 1=1,2, 



(4.13) 

Proof. Let us first consider the discontinuities of 8 1 (u(z) + A) and 6 2 (u(z) + A) across Ej. 
Let 7r : C — ► C be the projection of £ onto the Riemann sphere. Suppose z is a point in 
Ej. Let 2 G C and define the points z± e G Cj to be 



7r(4 ie ) = «±i C , z{ ie EC r (4.14) 

We will now choose e > to be real and positive and let z G Ej. From the definition of 
the integration contour in Figure [3] and the canonical basis of cycles in Figure 121 we have 
the following relation between the points z± ie as e — > 0. 



u\z 



"lie) 



U[Z 



Z G Z.g+1, 

9 



u(zl i£ ) = u(z 2 Ti£ ) + <f duj 



k=j 

9 



(4.15) 



u\z 



k=j 



J- 



j = 1, 



where e* is a vector with 1 in the k th entry and zero elsewhere and da; is the vector 

du = {du 1 ,...,duj g ) T . (4.16) 

From this and the periodicity of the theta function (14. 7B . we obtain 

6l(u(z) + A) = e 2 T (u(z) + A), zeEj, j = l,...,g + l. (4.17) 

Let us now consider a point z in Ej. Again, from the definition of the integration contour 
and the canonical basis, we have, as e — > 0, the following 



i+i 



u{z\ e ) = u(z l _ te ) + (-1) 

J \ i i i\M-i 



du, 



(4.18) 



//(:_) + (-iy +L Ue>, zeE v j 
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Figure 4: The contour E. 

where I = 1,2. From this and the periodicity of the theta function, we see that 

9 l + (u(z) + A) = 6 l _(u(z) + A)e ( " iy2 < uAz)+A ^) , zeEj, j = l,...,g. ( 4 - 19 ) 

From the definition of the integration contour, it is clear that 6 1 (u(z) + A) and 6 2 (u(z) + A) 
are analytic across R \ (Ai, \2g+2) and that 9 1 (u(z) + A) is analytic across S a -^ 2 . 

Let us now consider the discontinuities of 9 3 (u(z) + A) and 9 2 (u(z) + A) on S a ^^ 2 . Let 
z be a point on S , (7 _ A12 , from the definition of the contours, it follows immediately that 

6\{u{z) + A) = 9 3 _(u(z) + A), ze S a ., 2 . (4.20) 

Let us now consider the boundary value of 6 2 {u(z) + A) on the minus side of S a -^ 2 . For 
small and positive e — > 0, we have 



u(z 2 + e) + j> diUJ = u{z -e), (4.21) 

where £ is the close loop on C depicted in Figure HI Since this loop is contractible, we 
have 

6 2 _(u(z) + A) = 6 3 + (u(z) + A), ze S a ^ 2 . (4.22) 

Finally, the conditions 

el(u(z) + A) = 6 A T (u(z) + A), zeR, 

el( u (z) + A) = et(u(z) + A), zes,-,*. 

follow directly from the definition of the contour of integration. □ 
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4.2 Meromorphic differentials 

Another key ingredient in the construction of the parametrix is meromorphic differentials 
on the Riemann surface. Most of the results that we will be using can be found in [24]. 

Proposition 4. Let di, . . . , df. be k distinct points on a Riemann surface C. Let c±, . . . , Cf. 
be complex numbers with Y^=i c fc = 0. Then there exists a meromorphic 1-form dQ on C, 
holomorphic on C \ {di, . . . , d^} such that 

dQ(x) = + O(l)^ dxj, x^dj, j = l,...,k. (4.24) 

where Xj is a local coordinate near dj such that Xj = at dj . 

This result can be found for example, in [24] . (p. 52, theorem II. 5. 3) 
A meromorphic 1-form with simple poles only is called a meromorphic 1-form of the 
third type. Let dQ be a meromorphic 1-form of the third type. In order to define the 
periods of dfl unambiguously, one has to define the periods to be integrals around close 
loops dj and bj that are homo logic to the a and 6-cycles in Figure |2] in C \ fl po i e , where 
Q P oie is the set of poles of dfl 

By adding suitable multiples of holomorphic 1-forms to a given meromorphic 1-form, we 
can obtain meromorphic 1-forms with arbitrary a-periods. For example, if a meromorphic 
1-form given by Proposition H] has the following a-periods 



dQ = Aj, j = 1, 



Then the meromorphic 1-form 



dn = dn + J2(^j - A M 

will have a-periods 

dQ = Aj, j = l,...,g 



J a. 



but with the same pole structure and residues. Of course, we can not control both the a 
and the 6-periods of the 1-form. In fact, meromorphic 1-forms with prescribed a-period 
and pole structure is uniquely defined. A meromorphic 1-form with all a-periods zero is 
called a normalized meromorphic 1-form. 

Proposition 5. Let di, . . . , d& be k distinct points on a Riemann surface C. Let c±, . . . , Cf. 
be complex numbers with J^jLi c fc = and let Ax, . . . ,A g be arbitrary complex numbers. 
Then there exists a unique meromorphic 1-form of the third type dfl on C, holomorphic on 
£ \ {di, . . . , dfc} such that 

Res z=dj dQ = cj, j = l,...,k, 

a ■ i ( 4 - 25 ) 
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Proof. We have already shown the existence part. To see the uniqueness part, let dQ and 
dfi' be 2 meromorphic 1-forms of the third type with the properties (14.251) . Let dfl be their 
difference. Then, since both dQ and dtt' have the same singular behavior at the points 
dj, the 1-form dfl does not have any pole and is therefore holomorphic. Moreover, all its 
a-periods vanish. Since a holomorphic 1-form with vanishing a-periods has to be zero itself, 
(See, for example, [21], p. 65, Proposition III. 3. 3) the proposition is proven. □ 

We will conclude this section with a result that relates the periods of a normalized 
meromorphic 1-form to the values of the Abel map at its poles. 

Theorem 6. (See e.g. fEffl , p. 65, III. 3) Let rj be a meromorphic differential of the third 
type with simple poles at the points di G £ and fj be a holomorphic differential. Let IP and 
IP be their periods 



where xq is an arbitrary point on C 

5 Construction of the outer parametrix 

We will now construct the local parametrix with the the theta function and meromorphic 
1-forms. 

Let us now define a local coordinate w near oo 2 , the point on C 2 , £ 3 and £ 4 that 
projects onto oo in the Riemann sphere. 




(4.26) 



Then the Riemann bilinear relation is the following 




(4.27) 




OJZ 3 , 

2 -i 
OJ Z 3 



in the third quadrant of C 2 - 
in the first quadrant of £ 3 ; 




(5.1) 
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27TI , 1 



where u = e~ and the branch of z$ is chosen such that argz G (—11,11). One can check 
that w is indeed holomorphic in £ in a neighborhood of oo 2 . 

Let us now define four meromorphic 1-forms of the third type dAj, j = 1, . . . , 4 by the 
following properties. 

Definition 3. The normalized meromorphic 1-forms dAj are holomorphic in 

£ \ {±it, \{,..., Xl g+2} ti,..., i g , oo 1 , oo 2 }. 

where ±it are the points in £2 that project onto ±ic and ik are the zeros of 8(u(x)) . At 
these points they have simple poles with residues 

Res A idAj = Res-tiidAj = — -, k — 1, . . . , 2g + 2, 
k 1 

Res tfe dAj = 1, k = l,...,g, 

ReSooidAx = 0, Res oc 2dA 1 = 2, (5 2) 

Res oc idA 2 = 3, Res 00 2dA 2 = — 1, 

ReSooidAs = 2, Res oc 2dA 1 = 0, 

ReSooidA 4 = 1, ReSoo2dA 4 = 1, 

provided none of the i\ is equal to for some k. If some ti is equal to \\ for some k, then 
the residue at i\ will be ~. 

These 1-forms are then uniquely defined. We will denote the b-period of these 1-forms 
by /3j. 

/J 7 .= [/dA i ,...,/dA i ] , j = l,...,4. (5.3) 

\Jh Jbg J 

To avoid ambiguity in the b-periods, let it(ik) be the projection of ik on (See remark^. 
Then the b-periods are computed as integrals on b-cycles in £1 that intersects at any 
point x < ■n(ik) if vr(/ fc ) 7^ \ 2 k- If vr(i fc ) = X 2 k, then the b-cycle can intersect E k at any 
point x 7^ A 2 fc in C\. 

We will now define four functions in the Riemann surface £. First let 5^ e £ be the 
images of E k under the maps Ci,±( z ): that is, 



{(z,0\zeZ k , £ = £i, ± (s)}, k = l,...,g + l. (5.4) 



Let Zo be a point in 5 g+1 . The exact choice of zq is immaterial to the construction as long 
as z 7^ X\ g +i or X\ g +2- We will now define the functions Nj(z) on £ as follows. 

9 (u(z) + #^ + no) 
7VjUj " e 6(u(z)) (5-5) 

= e A ^ z) e j (z), a = ( ai ,...,a g ) T , j = l,...,4. 
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where the function Aj(z) is given by Aj(z) = f* dAj and the path of integration is defined 
in the same way as the ones for the Abel map, except that every path now starts at zq. 

Let z k be the point on Ck that projects to z in C. As before, we will now define four 
functions e A j ^ on the complex z-plane by 

Then these functions have the following jump discontinuities in the complex z-plane. 

Proposition 6. The functions e A ^ z ^ are analytic mC\(EU S a -^ 2 ). On K U S a -^, they 
satisfy the following conditions 



e A U^ 


= ^w, 


z G S fe , k — 1, 


...,0+1, 




= e Aj,_(z)+(- 




, fc = l,. 


e A i+(*) 


= e A U*\ 


z G \Eo U J , 


/ = 1,2, 


e A l±^ 


= ±e A i T W, 


Z G S a —^ 2 , 




e A l±^ 


= e A ,V*), 








= e A U*), 


2 G Sp—^ > I 1 


,4. 



/ = 1,2 



(5.7) 



where (f3j) k is the k component of the vector (3j. 

Proof. The proof follows similar argument as the ones used in the proof of Proposition [31 
First let us consider the jump discontinuities on S^. Let z be a point in and define the 
points z l ±ie as in f)4.14p in the proof of Proposition [31 First consider the boundary values 
e A J: +( z ) an d e A j-( z \ Choose integration contours X + and X_ from zq to the points z] t and 
z 2 ^ as in Figure El Let S = X + — E_. Then X can be deformed into the sum Ylf=k a i °f 
the a-cycles and a loop X 29+ 2 around the point ^l g +2 in £ \ A po i e , where A po i e is the set of 
poles of dAj. (See Figure E]). 

A po i e = {±it, A}, . . . , \l g+2 , Ll, ■ ■ ■ , Lg, OO 1 , OO 2 }. (5.8) 

Therefore we have 



exp (A,(4)) = exp ( A,-(z 2 ie ) + ^ / dA, + i" 



dAj . , 



(5.9) 



- exp (A^ziJ) , ^ G E fe , fc = 1, 



where the last equality follows from the fact that dA^ has residue — | at the point Xl g+2 - 

Let us now consider the boundary values e A j-^ and e A ^+^ on E&. Since zo G "g+i> 
the integration contours X_ and X + can now be chosen to lie in the lower (upper) half 
plane of Ci (£ 2 ). The loop X = X + — X_ can now be deformed into the sum — Yln=k a i °f 
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1 ic 


Ai A2 A3 A4 


<—ic \^ / 
S_ 



Figure 5: The contours £± for e A ^+ ( - 2 ' ) and e A ^-^. 
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the a-cycles. However, such a deformation will necessarily go pass the poles A^, . . . , X\ g+l 
and Lk, ■ ■ ■ ,i g of dAj (Recall that by the remark after Proposition [2} there is exactly one 
point Lk that belongs to E\ U Since the residues of dAj at these points are given by 
— 2 ( g ~ fc ) from the Xj and g — k from the ti when all ii and A^ are distinct, the total residue 
at these points is zero. It is clear from Definition [3] that, when some i\ coincide with the 
A^, the total residue at these points remains unchanged. Hence we have 

exp (A,(^J) = exp (A,(4) - jf dAj 

= exp (Aj(zD) , z G 2^, k=l,...,g. 

Let us now consider the boundary values on the gaps 5^. Let z G Efc. For the boundary 
values e A j>±^, we choose E± to be integration contours that go from zq to z± ie in £i 
without intersecting (— oo, X\ g+2 ) except at z and z± it . Let the loop S be S = S + — 
then for k = l,...,g, E can be deformed into the 6-cycle bk without passing any pole 
of dAj, except possibly Lk (Recall the definition of the 6-periods of dAj in Definition [3]). 
When Lk is not equal to X\ k or A2 fc+1 , dAj has integer residue at Lk, and when Lk is equal to 
either A^ or X\ k+l , the deformation from £ to bk will not have to go pass Lk- This implies 

e A}, +(2 ) = e Al_{ z )+^) kj ze z k , k = l,...,g. 

For k = 0, the loop S can be deformed into a small loop around the point oo 1 . Since the 
1-forms dAj have integer residues at oo 1 , we have 

If k = g + 1, then the loop S will be contractible in C \ A po [ e . Hence we have 

e ^, + W = e Aj,_W | zeE g+1 . 

On the other hand, for the boundary values e A ^ ±( - z ^ on H/-, let us consider S± to be inte- 
gration contours that go from zq to z± it in C 2 without intersecting (— oo, Xl g+2 ) except at 
z and z\ ifL . Let k = 1, . . . , g, and let vr(i fc ) be the projection of Lk onto the z-plane. Then 
depending on the relative positions of z, vr(t fc ) and A^, the loop S = S + — £_ can be de- 
formed into — 6^, together with small loops around the poles Lk, ■ ■ ■ ,L g and A2 fc+1 , . . . , X\ g+2 
in £ \ Ap i e ; or it can be deformed into the sum of — bk and small loops around the poles 
tfc+i, • • • ,L g and Xl k+1 , . . . , Xl g+2 in C \ A po i e . In either cases, the total residue of dAj at 
these points will be an integer. Therefore we have 



zeE k , k 



9- 



Similarly, for k — 0, the loop E can be deformed into a small loop around the points oo 2 
and ±it. Since the total residue the 1-form dAj at these points is an integer, we have 
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If k — g + 1, then the loop S will be contractible in £ \ A po i e . Hence we have 

Z E Eg +1 . 



We now consider the boundary values e A i>-^ and e A ^+^ at S^-^, let z E S a -^ 2 . Let us 
again denote by S + and XL contours of integration from zq to z — e in £3 and z + e in £2. 
Then depending on whether z is in the upper or lower half plane, the loop X = X + — XL 
can be deformed into to a small loop around the pole it or —it in £ \ A po i e (See Figure HI 
The loop X in this case is the same except that it begins and ends at z instead of X 2g+2 )- 



Since the residue of dAj around it or —it is — |, we have 



-2 G SL 



M2- 



(5.10) 



The rest of the jump discontinuities in ( 15. 7p now follow directly from the definition of the 
integration contours as in the proof of Proposition [3j □ 

Let us denote by iV- (z) the projection of Nj(z) onto the fc -sheet, that is, Nj(z) = 
Nj(z k ) = Nj {z,^k{z)), where £fc(z) is the function on Then we have the following. 

Theorem 7. Let N(z) be the 4x4 matrix whose elements are given by 



N jk {z) 



N"(z), 



Imz > 0; 



(-l) d ^Nf(z), lmz<0 ; 



(5.1i; 



where Nj(z) are defined in ( 15.51) . Suppose we have 



6 ( uioc 1 ) + ^ + na]e( u(oo 2 ) + A + net) ^ 0. (5.12) 



Let the constants Lj be 



Li 

L 3 



00 



Ln = lim N 



(W)W 



\rniN 3 -\w), 



lim 7V7 1 ' 



w }w, 



iu->0 



where w is the local coordinate near oo 2 defined in Ii5.1\) and the limits in L 2 . 
are taken as z — > oo 2 in the first quadrant of £ 2 ■ Then the matrix 



(I 








\ 





i 





IK 












i 





1° 











diag(L 1 ,L 2 ,L 3 ,L 4 ) N(z) 



(5.13) 
L 3 and L 4 

(5.14) 



satisfies the Riemann-Hilbert problem ( 3.7 ), where k is the following in the expansion of 



N 2 (z) at z 



00 



N 2 {z) 



L 2 1 (w 1 + K 2 fl 



KW 



as z — > oo 2 in the first quadrant of £ 2 . 
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Remark 3. The constants Lj, j = 1, 



Li 



, 4 can be represented as 



( lim e-^MuT 1 



#(«(oo 2 )) 



0( M (oo 2 ) + |^ + n«) 
#(u(oo 2 )) 



(5.15) 



lim e " A3(w) — 

Q (m(oo 2 ) + |^ + na) 



( lim e~ A4iw) w 



9(u(oo 2 )) 



9{u( 



A. 

27TJ 



+ no) 



where the limits are taken as z — > oo 2 m t/ie /irst quadrant of £2- 

Proof. First note that, by using Proposition [3] and Proposition EJ one can verify that N(z) 
does satisfy the jump discontinuities in (13.71) . 

Since N(z) satisfies the jump discontinuities of (13.71) . the matrix M(z)N~ 1 (z) does not 

2 

have any jump discontinuities in C. Moreover, this matrix does not grow faster than z$ 
at z = 00 and has at worst square root singularities at the points Aj and ±it. Since it 
has no jump discontinuities, all these singularities are removable and therefore we have 
M(z) = HN(z) for some constant matrix H. To determine the constant matrix H, we will 
have to study the behavior of N[z) as z — > 00. 
The behavior of M(z) is given by the following 



M{z) 



(l + 0{z- 1 ) O(z-I) o(z~ 2 r 
0{z- 1 ) * 
0(z- x ) * 
\ olz- 1 ) * 



* 
* 
* 



0(z~i 
* 
* 
* 



(5.16) 



where the 3x3 lower right block is given by 

/ -^(l + Oiz- 1 )) f^{l + 0{z- x )) 



^zHl + 0(z- x )) \ 



V3 



(i + o(z-i)) 



V3 



(l + 0(2Tf)) 



V3 

i 



\-^ 3 ZS (1 + 0(2-3)) ^-3(1 + 0(^3)) 



^-3(1+0(2-3 



(5.17) 



))/ 



for z — > 00 in the first quadrant. From the relation between the local coordinate w and z 
in ( 15.11) and the jump discontinuities of N(z) near oo 2 , we see that, if we can show that 
the functions Nj(z) behave as 



N l (z) = L- 1 \l + 0(z- 1 )), 
N 2 (z) = 0(z~ 1 ), z^oo 1 
N 3 (z) = Oiz' 1 ), z^oo^ 
N 4 (z) = O^" 1 ), z -> oo 1 



00 



N l {z) = 0{w 2 ), z 



00 



N 2 (z) = L2 1 (itf 

v 3 



KW + 0(w 2 )) , Z 



OO 



N 3 (z)=L 3 1 (l + 0(w 2 )), 



N 4 (z)=L^ L w(l + 0(w)), z 



00 



00 



(5.18) 
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when z — > oo 2 in the first quadrant of £2, then the matrix in (15.141) will be the unique 
solution of the Riemann-Hilbert problem (13.7p . The asymptotic behavior of Ni(z) and 
N±(z) follows immediately from the definition of the functions Nj(z) (15.51) . the constants 
Lj (I5.13P and behavior of the 1-forms dAj (15.21) . 

We will now prove the equations in (I5.18P for N 2 (z) and N 3 (z). Let the involution g on 
£ be g(z, £,(z)) = (—z,£(—z)). To simplify the notation, we shall simply denote g(z, £(z)) 
by —z. Let us consider the functions Nj(—z) for j = 2, 3. The singularity structure of this 
function is the same as Nj(z). By Proposition [3] and [6] and the expression of Nj(z) (15.51) . 
we see that the functions N 2 (z) and N 3 (z) satisfies the following jump discontinuities on 
£. 

N jt+ (z) = -Nj,_(z), zeE+, k = l,...,g + l, 

N jl+ (z) = e^ Mn ^N j ,-(z), zeE{, k = l,...,g, 1 = 1,2, (5.19) 



N h+ {z) = -N jt .(z), zeS: 



+ 



-H2 



where E k is defined in ( 15.41) . E l k is the interval on £[ that projects to E k . That is 

E l k = {(z,£)\zeE k , £ = 



The intervals S a _^ 2 are defined to be 

On the other hand, from (15.191) . we see that the function Nj(—z) has the following jump 
discontinuities 

N j>+ {-z) = -N^i-z), zeE~, k = l,...,g + l, 

N j>+ (-z) = e ( " 1)!+l2 ™^+ 1 - fe iV : ,-_(-z), zeE l k , k = l,...,g, 1 = 1,2, (5.20) 

Note that the union of the contours (u^H^) and (u^H^r) divides £ into 2 disjoint sets, 
which are the first sheet £± and the union of the other sheets £2 U £3 U £4. Similarly, the 
contour S^-^ U S^-m divides £ into the sets £1 U £ 2 and £ 3 U £ 4 . Let iV}^) be 

Nj{-z), z e £1 U £ 3 U £ 4 ; 



1 —NA—z), z e £ 2 . 



Since the constants a k satisfy the symmetry a k = 1 — (13.61) . from ( 15. 19ft . (15.20p and 

(15.211) we see that the function 

*,M - g| (5-22) 
is either a meromorphic function on £ with poles exactly at the g zeros of 

6 (u(z) + ^ + na 
\ 2m 
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(5.23) 



or it is a constant. By the assumption of the theorem, this theta function is not identically 
zero. Hence by Theorem El we see that Nj(z) must be a constant K.j. By using the 
jump discontinuities ( 15. 19ft . (15.20p of the Nj(z) near z = oo, and the relation between the 
coordinate z$ and w, we have the following behavior of Nj(z) near oo 2 

N 2 {z) = L 2 l (z* + k 2 , - kz~* + 0{w 2 ) 

N 3 (z) = L 3 1 (l + K 3fi z~^ + 0(w 2 )) , 

as z — > oo 2 in the first quadrant of C 2 , where the branch cut of z^ is chosen to be the 
negative real axis. On the other hand, since — z — ► oo 2 in the third quadrant when z — ► oo 2 
in the first quadrant, the functions Nj(z) have the following behavior 

iV 2 (z) = L 2 1 (-z* + K2)0 + KZ ~\ + 0(w 2 )) 
iV 3 (z) = L3 1 (l-^-l + O^ 2 ) 



(5.24) 

as 2 -> oo 2 in the first quadrant. On the other hand, since Nj(z) in (15.221) is a constant 
Kj, we also have 



N 2 {z) = JC 2 Lz l (z* + K 2fi - kz~^ + 0(w 
N 3 (z) = /C3L3 1 (l + k^z- 1 * + 0(w 



(5.25) 

as z — > cxd 2 in the first quadrant. By comparing (I5.24p and (I5.25p . we obtain (I5.18p . This 
concludes the proof of the theorem. □ 



6 The non-vanishing of the theta function 

We will now prove that the normalization constants 9 (^u(oo k ) + ^7 + > J ' = 1; • • • > 4 
and k = 1,2 does not vanish for any n G N. Then the solution S°°(z) of the Riemann- 
Hilbert problem (13.71) constructed in Theorem [7] exists and is well-defined. We will then 
show that it satisfies the conditions in Theorem HI We will use the results in Chapter 6 of 
[25]. 

First let us define a contour F that divides the Riemann surface £ into 2 halves. Let T 
be the set of points that is fixed under the map (f> in ( 14. 9p . That is, 

T = {xeC\ <p(x) = x} (6.1) 

Then T is a disjoint union of g + 1 closed curves Tj, j = 0, . . . , g on C, given by the 
followings. 

r = uf =0 r i5 

r = U 2 =1 (Eg U E k g+1 ) , (6.2) 
T j = U 2 k=1 E k j , j = l,...,g. 
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Figure 7: The dash lines indicate the loops Tj. 
where is the contour on £& that projects to Hj. That is 

s* = {(2,okG3 i) £ = &(*)}• 

In other words, the contours are the closed loops on £ that start from the branch point 
going through the interval [\2j,^2j+i] on C±, then enters C 2 at A^+i and go back to 
Xlj through the interval [\ 2 j, A27+1] on £ 2 . The contour T starts at X\, goes to —00 on 
the real axis on C%, then from +00 to Xl g+2 on the real axis on Ci, from which it enters 
£ 2 and goes to +00 along the real axis on £ 2 , then goes back from —00 on £ 2 to A \ along 
the real axis. (See Figure CJ). 

Note that the images of the cuts Sj on £1 and £ 2 do not belong to T. For example, let 
x = (z,£i i+ (z)) be a point on S], then 

<j>(x) = (z,£ 1& (z)) = (2, &,_(«)) ^ x. 

Similarly, the images of the real axis on £ 3 and £4 do not belong to T either. 

The curve T divides the Riemann surface C into 2 halves, C + and each of which 
is an open Riemann surface with boundary T. The Riemann surface C± consists of the 
upper (lower) half planes of £1, C 2 and £3 and the lower (upper) half plane of £4. The 
Riemann surface £ can now be thought of as a union of £ + , £_ and V. Moreover, the 
a-cycles defined in Figure [2] is homologic to the contours Tj. That is, we have 

r i~ a i> j = 1 ,---,9- (6-3) 
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We can think of C as the Riemann surface formed by gluing two copies of C + along the 
boundary V with an anti-holomorphic involution <f> that fixes V and maps £+ onto A 
Riemann surface formed in this way is called a Schottky double. Since £ is a Schottky 
double, we can apply the results in Chapter 6 of [25J to the theta function of C 
Let us define the tori § x and T x as in Propositions 6.2 and 6.8 of |25j . 

Definition 4. Let \ — (xi> • • • , Xg) T (Z/2Z) 9 and let Jo be the torus 

Jo = C 9 /A, A = Z 9 + Z 9 U. (6.4) 

The tori S x and T x are tori in Jo defined by 

S x = |sGJ ,| s = ^ X + ^, 

T x = {feJ ,l t = q + ~llx, ^R 9 

Note that this definition is different from the one in [25] because the theta function in 
[25] is defined differently. 

We can now apply the results in [25]. The first result tells us where the zeros Lj of the 
function 6(u(x)) are located. 

Proposition 7. (Proposition 6.4 of \EZ^) For any point Xq G T , s G S x , the function 
0(u(x) — u(xq) — s) either vanishes identically or has modulo 2, 1 + Xk zeros on Tk, where 
Xk is the k th component of the vector x- 

As a corollary, we have the following concerning the locations of the zeros ij. 

Corollary 1. The function 9(u(x)) has g zeros L\, . . . , i g such that i\. G Tk, k = 1, . . . , g . 

Proof. Let us take xq = Xl g+2 , X — and s = in Proposition [TJ then u(xq) = and 
by the paragraph after Lemma El we see that 8(u(x)) is not identically zero and hence by 
Proposition [71 it has 1 zero on each of the contour T^, k — 1, . . . ,g. □ 

The next result shows that the theta function does not vanish when its argument is 
real. 

Proposition 8. (Corollary 6.13 of \EM) Let T be the universal covering ofT , 

T = {f G C 9 , I F=?, ^i 9 .}. (6.6) 

Then the theta function 8(i) is real and positive for t £ To. That is, 6(i) is real and positive 
for allte R 9 . 

We can now prove that the periods j3j of dAj in Definition [3] are purely imaginary. 
This, together with Proposition [8] will imply the non-vanishing of the theta functions 

6 (u(oo k ) + §r + no) for k = 1, 2 and j = 1, . . . , 4. 
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^4 -^2g + 2 

1 1 — ic 



I 







1 ic 




*Ai '^^J^^u^ A* 


A 2g + 2 

i — ic 



Figure 8: The deformation from the cycle </>(ai) to a\ when ii G Sj. The deformations for 
other cycles are similar. 




Lemma 4. The periods f3j of the 1-forms dAj defined in Definition^ are purely imaginary. 

Proof. First note that, by Corollary [TJ all the points i\ and A; 1 are invariant under the 
involution <p. Hence the meromorphic 1-form dAj = dAj((f)(x)) has the same poles and 
residues as dAj(x). Let us show that the a-periods of dAj are zero. We have 

I dAj(0(x)) = I dAj(x) (6.7) 

From Figure [2J we see that the curve <f>(ak) consists of a path from the lower half plane 
in C\ that goes from Efc+i to and another path in the upper half plane of £2 that 
goes from S fc to S fc+1 . There are 3 poles of dAj between the loops a& and </>(a*.): A^ 
and A2 fc+1 (See Figure ED- From (15.21) . we see that dAj has a combined residue of at 
these points, and hence we can deform </>(a&) onto without affecting the value of (16.71) . 
Therefore, by (16. 7D . we see that 

<f dA,(0(x)) = O, j = 1, 

</a fc 
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By the uniqueness of normalized 1-form, this implies dAj(4>(x)) = dAj. Now we use ( 14. 10D 
for the 6-periods, since the relations for the 6-cycles in (" 14 .101) are exact and not up to 
deformation, we have 



dA^(x)) = -f dA,(x) = -(P 3 ) k . 

b k . Jb k 



where (13 j), is the k th component of the vector f3j. On the other hand, since dAj(<p(x)) = 
dAj, the above is also equal to ((3j) k . This implies the proposition. □ 

From Proposition [HJ and Lemma HI we obtain 

Theorem 8. There exists 5 > 0, independent on n, such that 9 (u(oo k ) + ^ + n&\ > 5, 
for k = 1, 2 and j — 1, ... ,4 and all n 6 N. 

Proof. Let us consider the normalized 1-form dQk that has simple poles at Xl g+2 and oo k 
with residues -1 and 1 for k — 1, 2. Then by similar argument used in the proof of Lemma 
HI we see that d£lk((j)(x)) = dflk and hence the 6-periods of dQk are all purely imaginary. 
Now by the Riemann bilinear formula (14.271) and the definition of the Abel map (14.81) . we 
see that 

2ttz (u l (oo k ) - ui(\\ g+2 j) = 2niu l (oo k ) = f dfi fc , k = 1,2, l = l,...,g. 

Jbi 

Hence u(oo k ), k = 1, 2 are real. 

Therefore, by Proposition [8] and Lemma HI we see that 6 (u(oo k ) + ^ + na) > 0. By 

the periodicity of the theta function (14. 7ft . we see that the theta function is in fact a map 
from T x R 9 -> C, where T is the torus T = R 9 /Z 9 . By Proposition [HI the restriction of 
the theta function on the compact set T x {0, 0, . . . , 0} is real and positive and hence there 
exists 5 > such that 9(t) > 5 for all t G T. This then implies the theorem. □ 

This implies that the function S°°(z) in Theorem [7] exists. We will now show that it 
satisfies the conditions in Theorem HI 

Corollary 2. The function S°°(z) in ^5.14 ) and its inverse (S* 00 ^)) -1 satisfy the conditions 



in Theorem^ 

Proof. Let us first show that the function N(z) in (I5.1ip is bounded in n uniformly in 
T, where T is defined in ( 13.91) . Since the entries of N(z) are restrictions of the functions 
Nj(z) in ( 15.51) on different sheets of the Riemann surface, we only need to show that Nj(z) 
is bounded inside the set T = Uf =1 ^(T) in C that projects onto T. From the periodicity 
property of the theta function ( 14. 7ft . we see that Nj(z) can be written as 

,(„(,) + & + na ) 

,W 6M) (6.8) 



A A Z ) i L 7 = 1 4 



6(u(z)) 
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where 7„ is a finite vector given by 

{in) i = na i ~ i na i}> 1 = !>•••> 9, 

where [x] is the biggest integer that is smaller than x. From (16.81) and the fact that 9(u(x)) 
is not identically zero (Proposition [5]), we see that Nj(z) is bounded in n uniformly in T. 

We will now show that the constants L 1; ...,L 4 in (15.151) are bounded in n. From 
the singularity behavior of the meromorphic 1-forms dAj in (15.21) . we see that following 
constants 

e- Al(ool) , lime"^^" 1 , lime- A3(,u) , lim e ~ A4{w)w , 

w— >0 w— >0 ui^O 

in (15.15P are all bounded and non-zero. Since they are all independent on n, they are also 
bounded away from infinity and zero as n — > oo. By Proposition [21 we see that 9(u(xj) is 
not identically zero and will only vanish at the points i\ that belong to IY Since neither 
oo 1 nor oo 2 belongs to T; for I = 1, . . . , g, the constants 6(u(oo k )), k — 1,2 are non-zero. 
Moreover, from the definition of the Abel map (14. §p . we see that ^(oo 1 ) and w(oo 2 ) are 
both finite and hence ^(^(oo 1 )) and 9(u(oo 2 )) are both bounded and are independent on 

n. Let us now consider the factors 9 (u(oo k ) + ^ + naj for k = 1, 2 and j = 1, . . . , 4. 

By Theorem [HI there exists 5 > 0, independent on n such that these constants are greater 
than 5. On the other hand, from the periodicity of the theta function (14 ,7p and the fact 
that the period matrix II is purely imaginary, (Lemma [3]) while the vector a in (13.51) is 

real, we see that 9 (^u(oo k ) + J^r + ndj is bounded in n as n — » oo. Hence the constants 

Li, . . . , L 4 in (15.141) and (15.151) are bounded away from infinity and zero as n — > oo. 

Finally, by considering the asymptotic expansion of Nj(z) in the local parameter w in 
(15. ip at z = oo and making use of (16.81) . we see that k in (15.141) is bounded in n as n — > oo. 
Since all the constants Lj and k are bounded in n as n — > oo and that all the Nj(z) are 
bounded in n uniformly in T, we see that S°°{z) is also bounded in n uniformly in T. To 
see that this is also the case for the inverse (S* 00 '(z))^ 1 , let us consider the determinant of 
S°°{z). Since S°°{z) is a solution to the Riemann-Hilbert problem (13.71) . the determinant 
det (S°°(z)) has no jump discontinuity in C and it behaves as 1 + 0(z~^) as z — > oo. From 
the expression of N(z) in (15. lip , we see that at Xj, only the first and second columns of 
N(z) have fourth-root singularities, while at the points ±zc, only the second and the third 
columns of N(z) have fourth- root singularities. Therefore the determinant of S°°(z) can 
at worst have square-root singularities at these points. Since det (S°°{z)) has no jump 
discontinuities in C, we see that det (S°°(z)) cannot have square- root singularities at these 
points. Hence det (S°°(z)) is holomorphic in the whole complex plane. By Liouville's 
theorem, this implies that det (S°°(z)) = 1. Since the entries of (S* 00 ^)) -1 are degree 3 
polynomials in the entries of S°°(z) divided by det (S°°(z)) = 1, we see that the entries of 
(S°°(z))~ 1 are also bounded in n uniformly in T. 

Finally, by considering the asymptotic expansion of Nj(z) in the local parameter w in 
(15. ip at z = oo and making use of (16.81) . it is easy to see that condition 2. in Theorem H] is 
satisfied for S°°(z) and its inverse. □ 
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We can now use Theorem H] to conclude that Theorem [2] and Theorem [3] are true. 
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